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The invariant yields, dN/dy, for J/ij) production at forward rapidity (1.2 < \y\ < 2.2) in U+U 
collisions at yJFffff=192> GeV have been measured as a function of collision centrality. The invariant 
yields and nuclear-modification factor Raa are presented and compared with those from Au-|-Au 
collisions in the same rapidity range. Additionally, the direct ratio of the invariant yields from U-l-U 
and Au-|-Au collisions within the same centrality class is presented, and used to investigate the role 
of cc coalescence. Two different parameterizations of the deformed Woods-Saxon distribution were 
used in Glauber calculations to determine the values of the number of nucleon-nucleon collisions in 
each centrality class, A^coii, and these were found to give significantly different Acoii values. Results 
using N„oii values from both deformed Woods-Saxon distributions are presented. The measured 
ratios show that the J/ij) suppression, relative to binary collision scaling, is similar in U-l-U and 
Au-|-Au for peripheral and midcentral collisions, but that J/'ij) show less suppression for the most 
central U-l-U collisions. The results are consistent with a picture in which, for central collisions, 
increase in the Jf'if yield due to cc coalescence becomes more important than the decrease in 
yield due to increased energy density. For midcentral collisions, the conclusions about the balance 
between cc coalescence and suppression depend on which deformed Woods-Saxon distribution is 
used to determine A^coii- 

PACS numbers: 25.75.Dw 


I. INTRODUCTION 

The study oi J/if production in high energy heavy ion 
collisions is motivated by the prediction that J/if forma¬ 
tion would be suppressed by color screening effects in the 
Quark Gluon Plasma (QGP) [T]. But relating J/if sup¬ 
pression to the energy densities of the hot matter formed 
in heavy ion collisions is complicated by the presence 
of competing effects that also modify J/if production. 
Gompeting effects [5] can be divided into cold nuclear 
matter (CNM) effects and hot matter effects. CNM ef¬ 
fects are those which modify the yield or kinematic dis¬ 
tributions oi J/if produced in a nuclear target in the ab¬ 
sence of a QGP. They include modification of the parton 
densities in a nucleus [Ml, breakup of the J/if or its cc 
precursor state in the nuclear target due to collisions with 
nucleons HHin], transverse momentum broadening due to 
the cc traversing the cold nucleus, and initial state parton 
energy loss m- CNM effects are expected to be strongly 
dependent on collision system, collision energy, rapidity 
and collision centrality. In hot matter, J/if yields can 
be enhanced by coalescence of cc pairs that are initially 
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unbound, but which become bound due to interactions 
with the medium |12j . At sufficiently high charm produc¬ 
tion rates, there can be a signihcant yield ol J/if from 
coalescence of a c and c from different hard processes. 

Precise J/if data extending down to zero pr have been 
published for Pb-|-Pb or Au-f An collisions at energies of 
= 17.3 GeV [B], 62.4 GeV [HUg, 200 GeV [Sl- 
[iHj and 2.76 TeV [T^. The nuclear-modification factor 
[Raa) of the J/if yield in the highest centrality colli¬ 
sions is observed to drop from = 17.3 GeV to 

200 GeV, and then increase strongly between 200 GeV 
and 2.76 TeV [T6| [ITl [19]. The rise in Raa between 
200 GeV and 2.76 TeV is well described [2^ in magni¬ 
tude and transverse momentum dependence by models 
that include coalescence of c and c pairs from different 
hard scattering processes [121 El]- The energy depen¬ 
dence of the modification suggests that J/if production, 
after accounting for the modihcation due to cold nuclear 
matter effects [lEdj, is increasingly suppressed from 
17.3 GeV to 200 GeV by stronger color screening in the 
increasingly hot QGP. But when the collision energy in¬ 
creases to 2.76 TeV, the rising underlying charm produc¬ 
tion rate leads to an increasing, and eventually dominant, 
contribution to the J/if cross section from coalescence. 

Of the collision energies observed so far the nuclear 
modification, Raa, for the most central collisions is at a 
minimum at y/s//// = 200 GeV, although there is little 
difference from 62.4 GeV, where the energy density is 
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smaller. The behavior of Raa in the range of energy 
densities accessed at energies between = 200 GeV 

in Au+Au collisions and 2.76 TeV in Pb+Pb collisions 
is not known. The measurement of J/ip yields in y/s~^ 
= 193 GeV 238u_|_238y collisions, the largest system yet 
studied at RHIC, during the 2012 data taking at the 
Relativistic Heavy Ion Gollider (RHIG) run provides an 
opportunity to increase the energy density above that for 
y'Sjyjy = 200 GeV Au+Au collisions by ^20% [52], and 
to observe the effect on the measured Raa ■ In this paper 
we report the results of measurements of the J/ij} yield 
in U+U collisions at forward and backward rapidity at 
y'SjyjY = 193 GeV using the PHENIX detector. 


II. PHENIX DETECTOR 


The U+U data used in this analysis were recorded at 
RHIC during 2012. The PHENIX detector [23] configu¬ 
ration used is shown in Fig. ^ The J/ip yields reported 
in this paper were all recorded in the muon spectrome¬ 
ters [5i] . 

The minimum-bias trigger was used for this data set. 
This required two or more hits in each arm of the beam 
beam counter (BBC), which comprises two quartz arrays 
of 64 Cerenkov counters, positioned symmetrically in the 
pseudorapidity range 3.0< \r]\< 3.9. The time difference 
between the two BBC arms provides the z-vertex posi¬ 
tion both in the minimum-bias trigger and in the offline 
analysis. 

The J/ip were reconstructed from their J/ij} -A 
decays, with the muons being detected in two muon spec¬ 
trometer arms that cover the rapidity ranges —2.2 < y < 
— 1.2 (south) and 1.2 < y < 2.4 (north). Each comprises 
a magnet, a copper and steel absorber followed by the 
muon tracker (MuTr) and the muon identifier (MuID). 
They are described in detail in Ref. [55]. An additional 
steel absorber material of thickness 36.2 cm was added in 
2010 to improve the muon yield relative to the hadronic 
background. The muon track momentum is measured in 
the MuTr, which has three stations, each comprising two 
or three chambers with two cathode strip planes each, 
with half of these planes having no stereo angle, and 
half having their cathode planes tilted at stereo angles 
that vary between 3 and 11.25 degrees. Discrimination 
between muons and punch-through hadrons is provided 
by the MuID, which comprises alternating steel absorber 
layers interleaved with larocci tubes. A muon candi¬ 
date is required to penetrate all the way to the last layer 
of the MuID, requiring a minimum muon momentum of 
3 GeV/c. The a cceptance for the T/V', discussed in detail 
in section mg is flat to within « 30% from transverse 
momentum, pT, of zero to 8 GeV/c. 



18.5 m= 60 ft 


FIG. 1. (Color online) A schematic side view of the PHENIX 
muon arms in 2012. The central arms, which are at midra¬ 
pidity, are not shown. 


III. DATA ANALYSIS 

For this analysis 1.08 x 10® events with primary ver¬ 
tices within ± 30 cm of the nominal interaction point were 
analyzed. This corresponds to an integrated luminosity 
of C =155.6 /ib“^ and a nucleon-nucleon integrated lu¬ 
minosity of / Cjyjydt = 238 x 238 x 155.6 ph~^ = 8.8 
pb“^. 

Because the invariant yields in the two muon arms 
must be identical for the symmetric U+U collision sys¬ 
tem, the analysis in the north muon arm was restricted 
to 1.2 < y < 2.2 to make the rapidity interval equal for 
the north and south arms. 


A. Centrality Determination 

The centrality determination for U+U collisions is 
based on the combined charge sum signals from the two 
BBC detectors. 

The U+U collisions were modeled using a Monte- 
Carlo simulation based on the Glauber model [55] with 
a deformed Woods-Saxon distribution for the U nucleus 
which accounts for its prolate nature. 


P = 


Po 

1 _|_ g([r-fl']/a) 


where R' depends on the polar angle 0: 


( 1 ) 


R' = i?[l + P2Y^{e) + p4Y!{9)\ (2) 

and where Y® is a Legendre polynomial, po is the normal 
nuclear density, R is the radius of the nucleus, and a is 
the surface diffuseness parameter. 

We considered two parameterizations of the deformed 
Woods-Saxon distribution for U. The first parameter set. 
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which we call set 1, is from Masui et al. m- The sec¬ 
ond, which we call set 2, is from Shou et al. [28]. The 
parameters for the two sets are summarized in Table [^ 
The parameterization of Shou et al. differs from the 
more conventional one of Masui et al. in two ways. First, 
the finite radius of the nucleon is taken into account. Sec¬ 
ond, rather than taking the mean radius and diffuseness 
for the deformed nucleus used in Equations and di¬ 
rectly from electron scattering experiments, their values 
are chosen so that after averaging over all orientations of 
the axis-of-symmetry for the nucleus, the average radius 
and diffuseness match the values reported from electron 
scattering experiments. The result is that the surface 
diffuseness, a, is considerably smaller for set 2 than for 
set 1, while the other parameters are similar in value. 


TABLE I. Parameters of the deformed Woods-Saxon distri¬ 
bution used in Equations and [^ 


Parameter 

set 1 

set 2 


R (fm) 

6.81 

6.86 


a (fm) 

0.6 

0.42 


02 

0.28 

0.265 


Si 

0.093 

0 



The smaller surface diffuseness of set 2 results in a 
notably more compact nucleus. The total inelastic U-fU 
cross section is 8.3 barns for set 1 and only 7.3 barns for 
set 2. The result is that the average number of binary 
collisions is 287 for set 1 and 323 for set 2. Although set 
2 appears to have a more consistent usage of the electron 
scattering data, we note that the neutron skin thickness 
of large nuclei is not probed via electron scattering and 
is thus rather unconstrained. 

The Glauber model was folded with a negative- 
binomial distribution, which models the charge produc¬ 
tion in the rapidity range of the BBC. The parameters 
of the negative-binomial distribution were fit to the mea¬ 
sured charge distribution from the BBC only in the signal 
range where the BBC minimum-bias trigger is known to 
be fully efficient. In the low signal range, the efficiency 
of the trigger was then determined from the ratio of the 
measured BBC charge distribution to the fitted-negative- 
binomial distribution. Figure [2 shows the distribution of 
measured BBC charge. Figure ^a) shows the measured 
charge distribution compared with the charge distribu¬ 
tion obtained from a Monte-Carlo calculation using the 
fitted-negative-binomial distribution. The efficiency of 
the BBC trigger was found to be 96±3%. Figure ib) 
shows the ratio of data to the Monte-Carlo calculation. 

Using the measured BBC charge, the events were di¬ 
vided into 10% wide centrality bins, as illustrated in 
Fig-i: a). Due to the limited statistical precision of the 
data sample, yields were extracted only for centralities 
from 0%-80%. The mean number of participating nu¬ 
cleons (IVpart) and mean number of nucleon-nucleon col¬ 
lisions (fVcoii) for each centrality bin were found from 



BBC Charge (North + South) 


FIG. 2. (Color online) The BBC charge distribution for 
U-l-U collisions for the north and south detectors combined, 
compared with a Monte-Carlo calculation using a negative 
binomial distribution that is fit to the data. The colored 
stripes show the BBC charge distribution divided into 10% 
wide centrality bins. 

the Monte-Carlo-Glauber calculation. In both cases the 
Glauber model used a nucleon-nucleon cross section of 
42 mb and assumed a hard core radius of 0.4 fm for the 
distribution of nucleons in the nucleus. The values of 
these parameters are summarized in Table [ll| for U-l-U, 
as determined using deformed Woods-Saxon parameter 
set 1 and 2. 

The systematic uncertainties for the mean A^part and 
fVcoii values in each centrality bin were estimated by vary¬ 
ing the parameters of the Glauber model within reason¬ 
able limits. The nucleon-nucleon cross section was varied 
from 42 mb down to 39 mb and up to 45 mb. For de¬ 
formed Woods-Saxon parameter set 1 the radius R and 
diffuseness parameter a were varied down to R = 6.50 fm 
and a = 0.594 fm and up to i? = 6.92 fm and a = 0.617 
fm. These variations were chosen to match the percent¬ 
age variations used for the Au-|-Au case, as were the vari¬ 
ations used in the evaluation of systematics for set 2. In 
addition to these variations of the Glauber model param¬ 
eters, the calculation was run without using a hard core 
for the nucleon distribution, and the trigger efficiency was 
varied within its uncertainty of 3%. For set 1, the un¬ 
certainty in the U deformation parameters has been esti¬ 
mated to be approximately 3% for the dipole component 
P 2 and approximately 50% for the small contribution of 
the quadrupole component /34 [55] • Varying j32 and /Si 
by these amounts resulted in an insignificant contribu¬ 
tion to the systematic uncertainty for Vpart and Vcoii. 
The systematic uncertainties are shown in Table |IT| 

The Vpart values obtained from the two deformed 
Woods-Saxon parameter sets are essentially indistin¬ 
guishable. However the differences between the N^oii val¬ 
ues obtained from the two parameter sets are, at some 
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TABLE II. Centrality parameters A^part and ALoii in U+U (this work) and Au+Au |16| collisions, estimated using the Glauber 
model. The systematic uncertainties are shown, and were estimated as described in the text. 



U-l-U set 1 

U±U set 2 

Au±Au 

Centrality 

Apart 

Acoll 

Apart 

Acoll 

A^part 

Acoll 

0%-10% 

386 ±5.2 

1161 ±126 

387 ±5.5 

1228 ±142 

326±3.9 

962 ± 97 

10%-20% 

273 ±6.7 

708 ± 69 

274 ±6.4 

770 ± 86 

236±5.5 

609± 60 

20%-30% 

190 ±6.6 

426 ± 42 

192 ±6.7 

473 ± 52 

168±5.8 

378± 37 

30%-40% 

128 ±6.7 

244 ± 24 

130 ±6.5 

277 ±29 

116±5.8 

224± 23 

40%-50% 

81.9±6.4 

130 ±16 

83.0 ±6.2 

149 ± 18 

76.2±5.5 

125± 15 

50%-60% 

47.7 ±5.4 

61.7±10 

48.5 ±5.3 

71.0± 11.3 

47.1±4.8 

63.9±9.4 

60%-70% 

24.7 ±4.0 

25.8 ±5.3 

25.3 ±3.9 

29.3 ±6.0 

26.7±3.7 

29.8±5.4 

70%-80% 

10.9 ±2.3 

9.2 ±2.4 

11.2±2.4 

10.2 ±2.8 

13.7±2.5 

12.6±2.8 


centralities, outside the uncertainties on the Ncoii values. 
This is due primarily to the large difference in the dif¬ 
fuseness values for the two sets. As noted earlier, the 
smaller diffuseness for set 2 (0.43 fm vs 0.6 fm for set 
1) combined with a similar mean radius results in larger 
A^coii values at all centralities because it makes the nu¬ 
cleus more compact. 

Because the A^coii values are important in the interpre¬ 
tation of the U-bU data, and are directly involved in the 
calculation of the Raa, we have chosen to consider the 
A^coii values from both sets 1 and 2 in the remainder of 
the paper. 


B. Muon-Track and Pair Reconstruction 

Muon candidates are charged particle tracks which 
penetrate all layers of the MuID. These MuID tracks are 
matched to MuTr tracks, which provide the momentum 
measurement. Although the requirement that tracks pass 
through the whole spectrometer arm reduces significantly 
the hadron contribution, a small percentage of charged 
hadrons can travel through without interacting (~0.1%), 
and these are a source of background. Some muons pro¬ 
duced from charged hadron decays in front of the MuTr 
are also reconstructed, and form a background of real 
muons in the spectrometer arms. Various offline analy¬ 
sis selection criteria are used to enhance the sample of 
good muon candidate tracks. There is a cut on the single 
track , and also on the difference in position and angle 
between the extrapolated MuTr and MuID parts of the 
candidate track. During the dimuon reconstruction, the 
selected track pair is fit with the event z-vertex position, 
and a second cut is applied for the track pair fit. 


C. Analysis 

The invariant mass distribution is formed by combin¬ 
ing all pairs of oppositely charged muon tracks. There 
is a significant combinatorial background under the J/^/> 


peak that is formed by single muons or mis-identified 
hadrons that randomly combine to form a pair. There 
is also a continuum due to correlated muons from semi- 
leptonic decays of open charm and bottom, and corre¬ 
lated muons from the Drell-Yan process. The combina¬ 
torial background may be estimated by event mixing, in 
which tracks with opposite charges from different events 
(of similar z-vertex and centrality) are combined ran¬ 
domly - see Ref. [T^ for more details. By construction, 
such mixing destroys any real muon correlations, and so 
the mixed event background does not reproduce the cor¬ 
related background from physics sources. 

Figure [^a)-(h) show the dimuon spectra from the 
south and north spectrometer arms for two exam¬ 
ple centrality bins and the combinatorial-background- 
subtracted mass spectra. The mixed-event combinatorial 
background is normalized to the real data using yields of 
muon pairs having the same charge sign found in both the 
mixed-pairs and real-pairs data samples in a range close 
to the J/if) mass peak region, 2.6< M< 3.6GeV/c^. The 
procedure is similar to that described in Ref. [3^. Be¬ 
cause the mass resolution precludes the separation of the 
J/ij} and peaks in this analysis, we present only J/ip 
results in this paper. 

The J/ij) yield was obtained using a fit in the mass 
range 1.7 to 5 GeV/c^. The fitting function included 
the mixed event background, a signal line shape, and 
an exponential. The signal line shape for the J/ip was 
a double Gaussian line shape modified by the mass de¬ 
pendence of the muon arm acceptance. The exponential 
was used to account for the correlated background that 
is not described by the mixed-event background, and to 
compensate for any systematic effects from under or over¬ 
subtracting the large combinatorial background. Because 
of the limited statistical sample for the U-|-U measure¬ 
ment, the form of the signal line shape was based on 
studies of the mass spectrum in p+p collisions, with in¬ 
put from studies performed for the Cu-|-Au system [31] 
where the J/ij: mass width was found to increase lin¬ 
early with increasing multiplicity in the muon arms. The 
functional form of the increase of the width with col- 
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FIG. 3. (Color online) Dimuon invariant mass spectra at forward and backward rapidity measured in the 0%-10% most central 
collisions (panels (a)-(d) )and in 60%-70% mid-peripheral collisions (panels (e)-(h)), integrated over the full pt range. Panels 
(a), (b), (e), and (f) show the invariant mass distribution, reconstructed from (solid symbols) same-event opposite charge-sign 
pairs and (open symbols) mixed-event pairs in U-l-U collisions. Panels (c), (d), (g), and (h) show the combinatorial background 
subtracted mass spectra. The solid line is a fit to the data using a double Gaussian line shape plus an exponential background, 
as described in the text. 


lision centrality was found independently for each muon 
arm from Cu-l-Au collision data, taken in the same RHIC 
data-taking period. This was then extrapolated where 
necessary to the higher multiplicity for the U-l-U sys¬ 
tem. As a cross check, a similar increase of the widths 
with multiplicity was observed in a GEANT [32] Monte- 
Carlo calculation with simulated pythia [33] J/ip events 
embedded into real U-l-U data events. The mass width 
varies from (peripheral to central) 0.14 to 0.2 GeV/c^ for 
the south arm, and from (peripheral to central) 0.14 to 
0.24 GeV !(? in the north arm. 

Because the width was fixed at each centrality by this 
procedure, the J/'ip yield was the only free parameter in 
the line shape. The fit is shown in Fig. He), (d), (g), 
and (h). The systematic uncertainty on the yield from 
the fit was estimated by varying the mass fit range, vary¬ 
ing the relative normalization of the signal to the mixed 
event combinatorial background by ±2%, and also ex¬ 
tracting the yield by using a like-sign combinatorial back¬ 


ground in the fit instead of the mixed-event background. 
Additional systematic checks were made by comparing 
the yields to those obtained from the raw signal count 
after exponential background subtraction, and to those 
obtained from other fit functions. The latter included al¬ 
lowing the Gaussian widths of the J/ip peak to be free in 
the fit function, where it was found that the yields agreed 
within the statistical uncertainty. The systematic uncer¬ 
tainty from the ht was larger in the north arm than in 
the south arm. For both arms, the fit systematic uncer¬ 
tainties were estimated to range from 0.8% for peripheral 
events to 10.9% for central events. 


D. Efficiency and Corrections 

The J/ij} detection efficiency was estimated in a two- 
step process. First, the efficiency to detect hits in each 
plane of the MuID was estimated by finding all roads 
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FIG. 4. (Color online) Acceptance (8) efficiency as a func¬ 
tion of collision centrality (Npart) for (squares) south arm, 
—2.2 < y < —1.2 and (circles) north arm, 1.2 < y < 2.2. The 
bands represent the uncertainty due the limited statistical 
precision of the embedding simulations. 


factors are reflected in the final data yields, where the sta¬ 
tistical uncertainties are largest for the north arm, and 
increase with increasing centrality. 


Between the measurement of the p+p data [34] used as 
the reference for the U-l-U nuclear modification (see sec¬ 
tion |iVA) and the U-|-U data, additional absorber ma¬ 
terial was added in front of the muon arms. The added 
absorber increases the minimum energy required for a 
muon to penetrate to the last gap of the MuID detector, 
and this reduces the J/'0 acceptance near y = 1.2 at low 
Pt- Because a realistic rapidity shape (from pythia) 
is used when calculating the rapidity integrated accep¬ 
tance 0 efficiency, it should be correct for both the U-l-U 
and p+p measurements. We note that the systematic 
uncertainty for the U-l-U acceptance 0 efficiency includes 
an uncertainty due to the possible deviation of the U-l-U 
rapidity shape from that given by pythia, as discussed 
in section IIII El 


made by charged particles through the MuID, but ig¬ 
noring hits in the plane of interest. The efficiency for 
the plane of interest was then estimated from the num¬ 
ber of roads with associated hits in that plane com¬ 
pared with all roads ignoring hits in that plane. Because 
the MuID efficiency decreases with increasing luminos¬ 
ity, the final efficiency used was the luminosity weighted 
average over all runs used in this analysis. The MuID 
efficiency is included in the calculation of the accep¬ 
tance 0 reconstruction efficiency described below. 

The full acceptance 0 efficiency (acceptance convo¬ 
luted with efficiency) is estimated by embedding pythia 
J/ip ^ pL~ decays - fully simulated via a GEANT de¬ 
scription of the PHENIX detector - into a reference data 
sample. The data plus simulation is then reconstructed 
with the same analysis as the real data and the hnal 
acceptance 0 efficiency is evaluated, normalized by the 
number of simulated J/ip ^ decays over the same 

rapidity range. 


TABLE III. Estimated systematic uncertainties. 


Source 

Uncertainty (%) 

Type 

J/ip signal extraction 

±0.8-10.9 

A 

input J/ijj Pt distributions 

±4.0 

B 

detector acceptance 

±5.0 

B 

reconstruction and trigger efficiency 

±5.0 

B 

run-to-run efficiency variation 

±2.8 

B 

Glauber (Vcoii) 

±10.8-33.0 

B 

p+p reference 

±7.1 

G 

p+p reference energy scale 

±3.6 

G 


The acceptance 0 efficiency for J/ij) in U-l-U collisions 
is shown in Fig. [^ As observed for other collision sys¬ 
tems, the north arm has lower efficiency than the south, 
and the efficiency is strongly centrality dependent. These 


E. Systematic Uncertainties 


The systematic uncertainties are divided into three 
groups: Type A - point-to-point uncorrelated uncertain¬ 
ties, Type B - point-to-point correlated uncertainties, 
and Type C - global scale uncertainties, which are sum¬ 
marized in Table jlllj The signal extraction systematic 
uncertainty due to the fitting procedure, discussed in 
Section |IIIC[ is treated as a Type A uncertainty. The 
uncertainty arising from the assumptions about the in¬ 
put J/tp momentum and rapidity distributions used in 
the PYTHIA calculations were previously studied in [14] 
and were found to be ^4.0%. The uncertainty in the de¬ 
tector acceptance is estimated to be ~5.0%. An overall 
uncertainty on the reconstruction and trigger efficiency 
obtained from embedding pythia events in real data is 
estimated to be ^5.0%. Small run-to-run variations in 
the MuID and MuTr efficiency were studied in m and 
found to be 2.8%. The uncertainty in the mean iVcou 
values for the centrality bins was studied as described in 
section III A and found to be 10.8-33.0% (see Table 0- 
The measured p+p reference invariant yield contributes a 
systematic uncertainty of 7.1%. This is smaller than the 
p+p global systematic uncertainty because both the BBC 
trigger efficiency for p+p and the iVcoii estimate used in 
the Raa depend on the assumed nucleon-nucleon cross 
section, and so their systematic uncertainties cancel in 
part when forming Raa- An additional small system¬ 
atic uncertainty is assigned to the p+p reference for the 
nuclear-modification factor because the p+p cross sec¬ 
tion, measured at ^/s = 200 GeV, had to be extrapolated 
to ^/s = 193 GeV to obtain the reference cross section 
for the U-l-U data set. That systematic uncertainty was 
taken conservatively to be equal to the correction, 3.6%, 
and when it is added in quadrature with the 7.1% from 
the p+p reference measurement, the overall global uncer¬ 
tainty due to the p+p reference becomes 8.1%. 
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IV. RESULTS 



FIG. 5. (Color online) Invariant yield of the J/'i/i at forward 
rapidity measured as a function of collision centrality (Vpart) 
for (closed circles) U+U, (open squares) Au+Au |16) . and 
(open diamonds) Cu+Cu |35| . 

The J/'0 yields measured in the two muon arms agree 
well, within statistical uncertainties. Therefore the U+U 
results presented here are averaged over both forward and 
backward rapidity. 


A. Yield and Raa versus Apart 


The measured invariant yield is obtained from Eq. 


B ^ 
dy 


I 


measured 


^event ^yAc 


(3) 


where B is the branching fraction for J/ij) decay to 

dimuons, .^measured measured number oi J/tp in¬ 

tegrated over all pT^ Ng^ent is the number of minimum- 
bias events analyzed, Ae is the acceptance 0 efficiency, 
and Ay is the rapidity range used when calculating the 
acceptance of the muon arms. 

The invariant yield for J/ip production in U+U colli¬ 
sions at ^/sNN = 193 GeV is shown versus Apart in Fig. 
where the vertical error bars represent the statistical un¬ 
certainty plus Type A systematic uncertainty for the 
yield extraction, added in quadrature. The boxes rep¬ 
resent the Type B systematic uncertainties, summed in 
quadrature. In addition to U+U, the largest system 
measured at RHIC, the Apart dependence of the invari¬ 
ant yield is shown for two other symmetric systems, 
Au+Au [1^ and Cu+Cu [35] at = 200 GeV. The 

yields below Apart ~ 150 are similar for all three systems 
at the same Apart- However for Apart ^ 200 the yield for 
U+U is larger than that for Au+Au. 

The departure from unity of the nuclear-modification 
factor. 


1 dN^^/dy 
~ (Aeon) dNPP/dy ’ 


(4) 


quantifies the modification of the Acoii normalized invari¬ 
ant yield in heavy ion collisions relative to the invariant 
yield in p-\-p collisions. The reference p-\-p invariant yield 
was obtained from [33] . Because the p+p data were mea¬ 
sured at y/s =200 GeV a scale factor of 0.964, determined 
from p+p PYTHIA simulations at ^/s =193 GeV and ^/s 
=200 GeV, was applied to the p+p invariant yield to ac¬ 
count for the difference in J/ip cross section. 
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FIG. 6. (Golor online) The nuclear-modification factor, 
Raa, measured as a function of collision centrality (Apart) 
for (closed circles) J/ip at forward rapidity in U+U collisions, 
for which the p+p reference data, measured at ,/s =200 GeV, 
has been scaled down by a factor 0.964 to account for the dif¬ 
ference in J/ip cross section between /~s =200 and 193 GeV. 
(open squares) Au+Au data [16]. (a) The U+U Raa calcu¬ 
lated using Woods-Saxon parameter set 1 m in the U+U 
Glauber model calculation, (b) The U+U Raa calculated 
using parameter set 2 |28 |. 

As discussed in section [HI A[ the deformed Woods- 
Saxon parameter sets 1 and 2 are derived using different 
assumptions, resulting in substantially different values of 
the surface diffuseness parameter, a. The authors of [21] 
argue that their approach in producing parameter set 2 
corrects deficiencies in the conventional method, and is 
the correct one. In any case, both assumptions cannot 
be correct, so the differences in the deformed Woods- 
Saxon parameters between the two sets do not represent 
a systematic uncertainty on the Raa, because one set 
is correct (within its uncertainties) and the other is not. 
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We have chosen to present the Raa calculated using both 
sets, so that the effect of using a conventional description 
and the description of reference [55] of the deformed U 
nucleus can be compared. 

The nuclear-modification factor for U-l-U collisions is 
shown as a function of A^part in Fig. § Fig.§ (a) shows 
the Raa calculated using A"coU values from deformed 
Woods-Saxon parameter set 1. Fig.[^(b) shows the Raa 
calculated using those from parameter set 2. The ver¬ 
tical bars represent the combined statistical and Type 
A systematic uncertainties and the boxes represent the 
Type B uncertainties. The Type C (global) uncertainties 
are listed in the legend. The overall global uncertainty 
(Type C) is 8.1%. The modification for Au-|-Au colli¬ 
sions [16] is shown for comparison. The Raa values for 
U-l-U collisions are provided in Table [TV] 

The measured Raa for the U-l-U collision system is 
similar to that for the Au-|-Au system, although the U-l-U 
data for the most central collisions show less modification 
than the Au-fAu data. 


B. Yield and Raa versus collision centrality 

Possible effects which may modify J/ip production in 
U-l-U collisions with respect to Au-|-Au collisions were 
discussed in |22|. The expected higher energy density 
in U-l-U compared to Au-|-Au (15%-20%) should lead to 
a stronger suppression due to color screening. On the 
other hand, for a given centrality, a larger N^oW value in 
U-l-U (compared to Au-|-Au) should increase charm pro¬ 
duction by cc coalescence. Cold nuclear matter effects 
due to shadowing are expected to be similar in both sys¬ 
tems. 



FIG. 7. (Color online) Invariant yield measured as a function 
of collision centrality for J/-)/) at forward rapidity for (closed 
circles) U-l-U and (open squares) Au-|-Au [16] collisions. 


As suggested in [55], we define for a given centrality 
class the relative nuclear-modification factor 


uu dN^^/dy f y 

AuAu dN^'^^'^/dy\0MixNgY\) ’ ^ 
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FIG. 8. (Color online) The relative nuclear-modification 
factor for U-l-U and Au-|-Au (Equation [^ as a function of 
collision centrality. Values were calculated using Aeon val¬ 
ues obtained from the U-l-U Glauber model calculation using 
Woods-Saxon parameter set 1 m and set 2 |28) . The values 
for set 1 and set 2 are slightly offset in centrality for clarity 
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where is the mean value of Ncoii in the given central¬ 
ity class (e.g. 10%-20%) in U-l-U collisions, and 
is the corresponding value for Au-|-Au collisions in the 
same centrality class. The N^^ii ratio in Equation is 
intended to account for the expected scaling of the J/ip 
cross section with in the case of cc coalescence. This 
assumes that the coalescence yield depends on the num¬ 
ber of charm quarks squared, and thus on the number 
of binary collisions squared. However to get the correct 
scaling in this case, the values of Ncoii for U-l-U have to 
be modified by a factor 0.964 because the charm produc¬ 
tion cross section in p+p collisions at =193 GeV is 
smaller than that at ^/s = 200 GeV by that factor. 

If the production of J/ip in central collisions was en¬ 
tirely due to cc coalescence, and if cold nuclear matter ef¬ 
fects were the same for both systems, the relative nuclear 
modification would be expected to be 1. The variable has 
the advantage that it is a direct ratio of U-l-U and Au-fAu 
invariant yields, eliminating the systematic uncertainties 
associated with the p+p reference when forming the Raa- 
The ratio of values for Au-|-Au and U-l-U also ap¬ 
pears in R-auAu- Fut the system atic uncertainty estima¬ 
tion for U-l-U (see section III AI was carried out in such a 
way that the systematic uncertainties in iVcou are highly 
correlated for U-l-U and Au-|-Au, and they mostly cancel 
when taking the ratio of Acoii values. There is still, how¬ 
ever, the large difference in the diffuseness parameter, a, 
coming from the difference between the models used to 
obtain parameter sets I and 2. We deal with that by 
calculating R^uAu both models. 


The invariant yield is shown plotted versus centrality 
class for U-l-U and Au-|-Au in Figure [7] To form the rel¬ 
ative nuclear modification of Equation]^ some rebinning 
of the Au-|-Au data is required. The rebinned invariant 
yields and their ratios are summarized in Table [^ Be¬ 
cause the U-l-U and Au-fAu data were taken in different 
years, and there were differences in the muon arm ab¬ 
sorber thickness due to an upgrade, all of the systematic 
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TABLE IV. The nuclear-modification factor (Raa) for U+U collisions (averaged over forward (1.2< y< 2.2) and backward 
(-2.2< y<-1.2) rapidity), as a function of centrality, derived from deformed Woods-Saxon parameter sets 1 and 2 (see text for 
details). The first and second uncertainties listed represent Type-A (statistical uncertainties plus point-to-point systematic 
uncertainties from the yield extraction) and Type-B uncertainties, respectively (see text for definitions). There is a Type-C 
(global) systematic uncertainty of 8.1%. 

Centrality 

Ruu (set 1) 

Ruu (set 2) 

0%-10% 

0.240 ± 0.045 ± 0.031 

0.225 ± 0.042 ± 0.031 

10%-20% 

0.314 ± 0.041 ± 0.035 

0.285 ± 0.038 ± 0.035 

20%-30% 

0.400 ± 0.028 ± 0.044 

0.349 ± 0.025 ± 0.044 

30%-40% 

0.463 ± 0.039 ± 0.053 

0.390 ± 0.033 ± 0.053 

40%-50% 

0.513 ± 0.040 ± 0.070 

0.410 ± 0.032 ± 0.070 

50%-60% 

0.604 ± 0.049 ± 0.109 

0.477 ± 0.039 ± 0.109 

60%-70% 

0.698 ± 0.077 ± 0.165 

0.567 ± 0.063 ± 0.165 

70%-80% 

0.646 ± 0.137 ± 0.241 

0.529 ± 0.104 ± 0.241 


uncertainties on the invariant yields are propagated in 
the ratio. 



FIG. 9. (Color online) The ratio of the invariant yield for 
U+U to that for Au+Au, as a function of collision centrality. 
The curves show how the ratio would vary with centrality if 
the invariant yield scaled with Acoii (dashed curve) or with 
A^coii^ (solid curve). The A^coii ratio curves obtained from the 
U+U Glauber model calculation using Woods-Saxon param¬ 
eter set 1 m are shown in blue, those using parameter set 
2 are shown in red. The Acoii value for U+U is multi¬ 
plied by a factor of 0.964 to account for the decrease of the 
charm cross section for p+p collisions from 200 GeV to 193 
GeV collision energy. 


on Nco\i- For comparison, in Figure we also present 
curves showing how the ratio would depend on centrality 
if it scaled with iVcoii (dashed curve) or with Ncoii^ (solid 
curve). As in Equation]^ the values of Acoii for U+U are 
multiplied by 0.964 to account for the difference in cross 
section for J/'tj) production in p+p collisions at 193 GeV 
and 200 GeV collision energy. When Ncoii values from set 
2 are used, the measured ratios for midcentral collisions 
lie at or below the ratio of iVcoii values, but as the collision 
centrality increases the data points move above the ratio 
of Vcoii values until, for the most central collisions, they 
favor the Vcoii^ curve. When Ncoii values from set 1 are 
used, the data slightly favor the trend of the Vcoii^ curve 
across the centrality range. 

These comparisons of the ratios of the U+U and 
Au+Au invariant yields with Vcoii values derived from 
both deformed Woods-Saxon parameter sets suggests 
that, for central collisions, increased cc coalescence is 
more important than stronger suppression due to in¬ 
creased energy density. When the Vcoii values from set 2 
are used, the peripheral and midcentral data suggest that 
the increased suppression due to increased energy density 
and increased coalescence due to larger underlying charm 
yields approximately cancel each other. When the Ncow 
values from set 1 are used, the peripheral and midcentral 
ratios are consistent with either N^oii or Vcou^ scaling. 


The values of R^^Au plotted in Figure]^ and sum¬ 
marized in Table |vg For the N^ow values from set 2, the 
relative nuclear modification falls below one for collisions 
in the 40%-60% region, but rises to one for the most cen¬ 
tral collisions. For Vcoii values from set 1, the relative 
nuclear modihcation is slightly above 1 for central colli¬ 
sions and approximately one across the remainder of the 
centrality range. These results suggest that the invariant 
yield scales with Vcoii^ for the most central collisions, at 
least. 

As a different way of looking at the data, we show 
in Figure the ratio of the invariant yields for U+U 
and Au+Au, taken from Table [V] These do not depend 


V. SUMMARY 

We have presented measurements of the pp integrated 
invariant yield dN/dy and nuclear modification Raa for 
J/ijj production in U+U collisions at ^Jsnn = 193 GeV, 
and compared them with existing data for Au+Au col¬ 
lisions at -xJsmn — 200 GeV. In addition to comparing 
the invariant yields and nuclear modification, Raa, for 
the two systems, we have combined them to form the 
relative nuclear modification, .RaJJau 122] ■ The relative 
nuclear modihcation (Equation!^ was proposed as a way 
to eliminate the systematic uncertainties associated with 
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TABLE V. The invariant yield for U+U collisions (averaged over forward (1.2< y< 2.2) and backward (-2.2< y<-1.2) rapidity), 
the invariant yield for Au+Au collisions, and the ratio of invariant yields plotted in Fig. all as a function of centrality. The 
first and second uncertainties listed represent Type-A (statistical uncertainties plus point-to-point systematic uncertainties from 
the yield extraction) and Type-B uncertainties, respectively (see text for definitions). There is no Type-C (global) systematic 
uncertainty for the invariant yield. 


Centrality 

B /dy 


/dy 

0%-10% 

182.7+ 35.1 ± 15.7 

105.3 ± 18.8 ±1^4 

1.73 ± 0.46 ±E!:i 

10%-20% 

144.5 ± 19.8 ± 12.4 

92.3 ± 10.7 ±1111:^ 

1.57 ± 0.28 ± 0.22 

20%-30% 

110.6 ± 8.2 ± 9.5 

75.8 ± 6.3 ± 6.4 

1.46 ± 0.16 ± 0.18 

30%-40% 

72.4 ± 6.3 ± 6.2 

56.8 ± 3.8 ± 5.1 

1.27 ± 0.14 ± 0.16 

40%-50% 

41.7 ± 3.4 ± 3.6 

40.5 ± 2.4 ± 3.4 

1.03 ± 0.10 ± 0.12 

50%-60% 

23.8 ± 2.0 ± 2.0 

28.2 ± 1.4 ± 2.1 

0.84 ± 0.08 ± 0.10 

60%-70% 

11.9 ± 1.4 ± 1.0 

13.4 ± 0.8 ± 1.0 

0.89 ± 0.12 ± 0.10 

70%-80% 

4.0 ± 0.8 ± 0.4 

7.9 ± 0.54 ± 0.60 

0.51 ± 0.11 ± 0.06 


TABLE VI. The relative nuclear-modification factor for U-l-U 
and Au+Au collisions (see Equation]^ (averaged over forward 
(1.2< y< 2.2) and backward (-2.2< y<-l.2) rapidity). The 
values derived from the deformed Woods-Saxon parameter 
sets 1 [23 and 2 |28| are shown. The first and second uncer¬ 
tainties listed represent Type-A (statistical uncertainties plus 
point-to-point systematic uncertainties from the yield extrac¬ 
tion) and Type-B uncertainties, respectively (see text for defi¬ 
nitions). There is no Type-C (global) systematic uncertainty. 


Centrality ^AuAu (set 1) RauAu (set 2) 


1 

0 

10 % 

1.29 

1 i 

-H 

CO 

0 

0.23 

:0.16 

1.13 

1 ± 

: 0 . 3 C 

1 i 

0.20 

:o.l4 

10 %- 

- 20 % 

1.27 

± 

0.23 

± 

0.18 

1.04 

± 

0.19 

± 

0.15 

20 %- 

- 30 % 

1.26 

± 

0.14 

± 

0.15 

0.96 

± 

0.11 

± 

0.12 

30 %- 

-40% 

1.21 

± 

0.13 

± 

0.15 

0.86 

± 

0.09 

± 

0.11 

40 %- 

- 50 % 

1.12 

± 

0.11 

± 

0.13 

0.72 

± 

0.07 

± 

0.09 

50 %- 

- 60 % 

1.03 

± 

0.10 

± 

0.12 

0.64 

± 

0.06 

± 

0.07 

60 %- 

- 70 % 

1.25 

± 

0.16 

± 

0.14 

0.83 

± 

0.11 

± 

0.10 

70 %- 

- 80 % 

0.96 

± 

0.21 

± 

0.11 

0.64 

± 

0.14 

± 

0.07 


the formation of Raa, and to cancel most of the system¬ 
atic uncertainties associated with the number of binary 
nucleon collisions, Ncoii- It is designed to have a value of 
one if the J/ip cross section is dominated by cc coales¬ 
cence. We have also compared the ratio of the invariant 
yields for U+U and Au+Au with the ratio of A^coii values 
and the ratio of Acon^ values for U+U and Au+Au. 

We discussed the effect of using two different param- 
eterizations for the deformed Woods-Saxon distribution 
on the estimates of A^coii for U+U collisions. A recently 
proposed method for estimating the Glauber model pa¬ 
rameters (set 2) [3S] leads to a smaller surface diffuseness 
for U, and thus larger values of iVcoii, than a conventional 
estimate (set 1) [57]. We presented Raa values and rela¬ 
tive nuclear modification values obtained using N^oW val¬ 
ues from both deformed Woods-Saxon parameter sets. 

For both sets of N^ow values the Raa for U+U is found 
to be less suppressed than for Au+Au in central colli¬ 


sions that have a similar number of participants. The 
relative nuclear modification is found to be one for the 
most central collisions for both NcoW sets, but for set 2 
it falls below one for mid-peripheral and peripheral colli¬ 
sions. When the ratios of invariant yields are compared 
with the ratios of N^ow and A^coii^ values, they are found 
to show a slight preference for Ncovi^ scaling for central 
collisions. 

For both sets of NcoW values the behavior is consis¬ 
tent with a picture in which, for central collisions, the 
increase in J/ij; yield for U+U due to cc coalescence be¬ 
comes more important than the decrease in yield due to 
increased energy density. For WoU values from set 1, the 
results are consistent with both A^coii and N^ow^ scaling 
for mid-peripheral collisions. For iVcoii values from set 2, 
the results suggest that in the 40%-60% centrality range 
the increased suppression due to higher energy density in 
U+U collisions is more important than the increased J/^/> 
yield due to coalescence caused by the higher underlying 
charm production. 
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